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1.  Introduction 


In  Contract  DAAD19-99-C-004,  numerical  simulation  of  projectile  base  flows  yielded  promising  results. 
Differences  between  computed  and  measured  base  pressure  and  reattachment  points  over  a  range  of  Mach  and 
Reynolds  numbers  were  acceptably  small  for  most  design  studies  [Wilcox  (2001)].  On  balance,  a  modified 
version  of  the  Wilcox  Stress-co  model  [Wilcox  (1998)]  predicted  averaged  base  pressure  a  bit  closer  to 
measurements  than  the  k-co  model. 


However,  in  order  to  achieve  the  improved  results  with  the  Stress-co  model  and  to  improve  numerical  stability, 
we  modified  the  production  term  in  the  co  equation.  Specifically,  we  made  the  following  replacement. 


CD  dut 

aTr'*ree's' 


In  Equation  (1),  a  is  a  closure  coefficient,  co  is  specific  dissipation  rate,  ry  is  the  specific  Reynolds-stress  tensor, 
Ui  is  the  velocity  vector,  xj  is  the  position  vector  and  Sy  is  the  mean  strain-rate  tensor.  We  made  no  assessment 
of  the  proper  value  for  the  closure  coefficient  a  in  Contract  DAAD19-99-C-004. 

This  study  accomplished  five  things. 

1.  It  focused  on  the  appropriateness  of  the  model  modification  quoted  in  Equation  (1). 

2.  It  examined  details  of  one  of  the  computations  done  in  Contract  DAAD19-99-C-004. 

3.  It  did  a  grid-resolution  study  to  confirm  the  accuracy  of  the  computation  analyzed  in  detail. 

4.  It  attempted  improvements  to  the  Stress-co  model  and  tested  a  nonlinear  stress/strain-rate  constitutive 
relationship  with  the  k-co  model. 

5.  It  made  computations  for  effects  of  base  bleed. 

Section  2  describes  the  first  part  of  our  analysis  addressing  the  question  of  the  appropriate  value  for  a.  We 
discovered  that  the  surface  boundary  condition  for  both  the  Stress-co  and  k-co  models  can  be  improved  in  a  way 
that  improves  the  accuracy  of  both  models  for  boundary  layers.  We  also  discovered  that  no  choice  of  the 
coefficient  a  exists  that  warrants  the  change  in  the  production  term. 

Section  3  includes  detailed  plots  of  velocity  and  turbulence-property  profiles  and  comparison  with 
corresponding  measurements  for  the  Herrin-Dutton,  Mach  2.5  base  flow  experiment  [Herrin-Dutton  (1994)].  It 
also  summarizes  results  of  a  grid-resolution  study. 

Section  4  discusses  attempts  to  improve  the  Stress-co  model  and  to  implement  a  nonlinear  stress/strain-rate 
constitutive  relationship  with  the  k-co  model. 

Section  5  presents  new  results  in  which  base  bleed  has  been  included. 


2.  Surface  Boundary  Conditions 

Whenever  a  turbulence  model  is  modified,  it  should  be  retested  for  turbulent  flow  data  against  which  the  model 
was  originally  calibrated.  One  set  of  data  for  which  both  the  Stress-cu  and  k-co  models  have  been  calibrated  is 
the  viscous  sublayer.  In  particular,  in  analyzing  the  sublayer,  we  have  devised  surface  boundary  conditions  for 
both  models  that  are  appropriate  for  both  smooth  and  rough  surfaces. 

As  noted  in  the  Introduction,  we  found  in  Contract  DAAD19-99-C-004  that  the  Stress-cn  model’s  predictions 
and  numerical  stability  for  base  flows  improve  by  revising  the  production  term  in  the  co  equation.  However, 
while  the  production-term  revision  appeared  to  be  an  improvement  to  the  model,  we  had  not  verified  that  the 
revision  is  satisfactory  for  all  of  the  other  flows  to  which  the  Stress-cu  model  has  been  applied. 

In  this  research  project,  we  found  from  study  of  the  viscous  sublayer  and  attached  boundary  layers  that  the 
closure  coefficients  would  have  to  be  a  bit  different  from  those  used  in  the  original  model.  In  doing  the 
analysis,  we  made  a  remarkable  discovery  regarding  surface  boundary  conditions. 

The  original  Stress-cu  model  used  surface  boundary  conditions  that  are  somewhat  different  from  those 
appropriate  for  the  k-co  model  [Wilcox  (1998)].  In  analyzing  the  new  version  of  the  Stress-cu  model,  we 
conclude  that  both  models  should  use  exactly  the  same  boundary  condition  for  both  smooth  and  rough  surfaces. 
Specifically,  for  both  models,  the  surface  value  of  co  is  given  by 


at  y  -  0 


(2) 


where  the  dimensionless  quantity  Sr  is  the  following  function  of  the  dimensionless  roughness  height  defined  by 
ks+  =  uTks/v  ( ks  is  the  physical  surface-roughness  height). 
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For  very  smooth  surfaces,  we  can  deduce  from  these  equations  that  the  appropriate  surface  boundary  condition 
(the  so-called  “slightly-rough-wall”  boundary  condition)  is  given  by 
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The  figure  below  shows  the  effect  of  surface  roughness  on  model-predicted  skin  friction,  c/,  for  boundary  layers 
at  both  Mach  5  and  for  incompressible  conditions.  As  shown,  consistent  with  measurements,  there  is  no 
noticeable  effect  of  roughness  when  ks+  is  less  than  5.  The  results,  which  are  virtually  identical  for  the  k-co  and 
Stress-cu  models,  are  clearly  consistent  with  the  measurements  of  Reda,  Ketter  and  Fan  (1974). 
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Figure  1:  Computed  and  measured  effects  of  surface  roughness  on  skin  friction  for  compressible  turbulent 
boundary  layers. 


Interestingly,  the  boundary  condition  described  in  Equations  (2)  -  (4)  is  essentially  independent  of  the  revised 
production  term  in  the  co  equation. 

Moving  on  to  attached  boundary  layers  and  homogeneous  turbulent  flows,  it  quickly  became  obvious  that 
changing  the  production  term  as  described  in  Equation  (1)  had  an  undesirable  effect  on  model  predictions  for 
these  baseline  test  cases.  No  value  of  the  closure  coefficient  a  exists  that  preserves  the  model’s  integrity. 
Consequently,  the  modification  described  in  Equation  (1)  must  be  rejected. 


3.  Detailed  Analysis  of  a  Mach  2.5 
Projectile  Base  Flow 

A  central  task  of  this  project  was  to  analyze  results  of  the  k-co  and  Stress-cu  model  computations  performed  for 
base  flows  in  Contract  DAAD19-99-C-0047.  Specifically,  the  goal  has  been  to  make  plots  throughout  the  near¬ 
wake  region  of  the  Herrin-Dutton,  Mach  2.5  base  flow  and  compare  computed  and  measured  flow  properties. 

Figure  2  compares  computed  and  measured  horizontal-velocity  profiles,  u/U^,  at  four  stations  downstream  of 
the  base.  Table  1  summarizes  the  measurement  locations  and  the  computed  and  measured  centerline  stagnation 
point  locations. 


Table  1.  Measurement  Stations  and  Stagnation-Point  Locations 


Item 

x/D 

Measurement  Station  1 

1.26 

Stagnation  Point:  k-co 

1.29 

Stagnation  Point:  Measured 

1.32 

Measurement  Station  2 

1.42 

Stagnation  Point:  Stress-ry 

1.43 

Measurement  Station  3 

1.73 

Measurement  Station  4 

1.89 

Stagnation  Point:  RNG  k-e 

1.93 

3.1  Horizontal- Velocity  Profiles 

As  shown  in  Figure  1,  the  k-co  model’s  horizontal  velocity  profiles  are  quite  close  to  measured  values  at  the  first 
two  Measurement  Stations.  This  is  consistent  with  the  fact  that  these  two  profiles  straddle  the  measured 
centerline  stagnation  point,  and  that  the  k-co  model’s  stagnation  point  also  lies  between  Measurement  Stations  1 
and  2.  By  contrast,  the  Stress-m  model’s  stagnation  point  lies  just  downstream  of  Measurement  Station  2  so 
that  greater  discrepancies  are  expected. 

Both  models’  predictions  are  far  closer  to  measurements  than  the  RNG  k-e  model  as  implemented  by  Sahu 
(1994).  The  RNG  k-e  model  predicts  a  stagnation  point  that  lies  downstream  of  all  four  Measurement  Stations, 
corresponding  to  a  recirculation  region  much  greater  in  extent  than  measured. 

At  the  two  Measurement  Stations  farther  downstream  of  the  base,  the  k-co  model’s  predicted  horizontal 
velocities  are  1 5%-20%  smaller  than  measured,  while  those  of  the  Stress-co  model  are  20%  to  25%  smaller  than 
measured.  This  is  typical  of  RANS  turbulence-model  predictions  for  reattaching  flows,  i.e.,  the  computed 
approach  to  a  new  equilibrium  state  is  slower  than  observed  experimentally. 
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Figure  1.  Comparison  of  computed  and  measured  horizontal  velocities  downstream  of  the  base  for  the 
Herrin-Dutton,  Mach  2.5  flow  past  an  axisymmetric  body  with  a  square-cut  base. 


3.2  Radial  Velocity  Profiles 

Figure  2  compares  computed  and  measured  radial-velocity  profiles.  The  most  noteworthy  difference  between 
computed  and  measured  flowfields  is  shown  in  the  location  of  the  reattachment  shock.  The  measured  profiles  at 
x/D  =  1.26  and  x/D  =  1.42  both  show  a  sharp  change  in  the  magnitude  of  the  radial  velocity  near  r/R  =  0.6. 
Farther  downstream  at  r/R  =  1 .73  the  rapid  change  occurs  near  r/R  =  0.9,  while  it  is  off  scale  at  x/D  =  1.89. 


By  contrast,  the  k-co  model  fails  to  predict  abrupt  changes  in  radial  velocity.  Although  the  Stress-m  model  does 
predict  sudden  changes  in  radial  velocity,  the  changes  occur  at  about  twice  the  measured  distance  from  the 
centerline.  Since  the  S  tress-tv  model  computation  required  finer  resolution  near  the  centerline  than  the  k-co 
model  computation,  the  absence  of  abrupt  changes  in  radial  velocities  suggests  the  possibility  of  inadequate 
resolution  of  the  flowfield.  Inspection  of  contour  plots  clearly  indicates  formation  of  a  reattachment  shock  for 
the  k-co  model  computation. 
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Figure  2.  Comparison  of  computed  and  measured  radial  velocities  downstream  of  the  base  for  the  Herrin- 
Dutton,  Mach  2.5  flow  past  an  axisymmetric  body  with  a  square-cut  base. 


3.3  Turbulence  Profiles 


Figures  3  and  4  compare  computed  and  measured  turbulence  kinetic  energy,  k/uJ,  and  Reynolds  shear  stress, 

2 

<-u  V  ’>/  UJ,  throughout  the  near-wake.  Both  models  fail  to  predict  turbulence  kinetic  energy  profiles  that  are 
close  to  the  measured  profiles.  However,  inspection  of  the  Reynolds  shear  stress  profiles  reveals  a  remarkable 
result.  The  k-co  model’s  shear  stresses  are  much  closer  to  measured  stresses  than  those  of  the  Stress-co  model. 
This  is  the  reason  the  velocities  are  closer  as  well. 

As  will  be  discussed  in  Section  4  below,  rerunning  the  Stress-co  model  computations  with  several  adjustments 
in  closure  coefficient  values  and  surface  boundary  conditions  yields  no  significant  reduction  in  differences 
between  computed  and  measured  flow  properties. 

Nevertheless,  both  models  yield  a  large  improvement  in  predictive  accuracy  relative  to  the  k-s  model.  While 
base  pressure  varies  with  radius,  the  average  value  is  within  a  few  percent  of  the  measured  value.  The  model 
thus  provides  an  accurate  prediction  of  body  drag.  Also,  the  reattachment  length  differs  from  the  measured 
length  by  just  2%. 

3.4  Grid-Resolution  Study 

To  investigate  the  possibility  that  the  k-co  computation  might  have  insufficient  grid  resolution,  we  repeated  the 
computation  with  four  times  as  many  grid  points.  The  original  computation  used  200  grid  points  in  the 
streamwise  direction  and  160  points  between  lower  and  upper  mesh  boundaries.  We  generated  a  grid  with  400 
streamwise  points  and  320  points  in  the  other  direction. 

Changes  in  major  structural  features  such  as  base  pressure  and  location  of  the  centerline  stagnation  point  are 
hardly  distinguishable.  Richardson  extrapolation  indicates  that,  overall,  the  numerical  results  are  within  3-4% 
of  the  continuum  solution.  Hence,  we  conclude  that  the  results  shown  in  Figures  1-4  are  accurate. 
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Figure  3.  Comparison  of  computed  and  measured  turbulence  kinetic  energy  profiles  downstream  of  the  base 
for  the  Herrin-Dutton ,  Mach  2.5  flow  past  an  axisymmetric  body  with  a  square-cut  base. 


Reynolds  Shear  Stress  at  x/D  =  1.42 
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Figure  4.  Comparison  of  computed  and  measured  turbulence  Reynolds  shear  stress  profiles  downstream  of 
the  base  for  the  Herrin-Dutton ,  Mach  2.5  flow  past  an  axisymmetric  body  with  a  square-cut  base. 


4.  Turbulence-Model  Modifications 


In  an  attempt  to  improve  model  predictions,  we  tried  two  things.  First,  for  the  Stress-tv  model,  we  did 
numerical  experiments  with  one  of  the  model’s  key  closure  coefficients.  Second,  we  used  a  nonlinear 
Reynolds-stress/strain-rate  constitutive  relation  with  the  k-co  model.  Neither  modification  proved  fruitful. 

4.1  Stress-co  Model  Revision 

One  of  the  terms  in  the  equation  governing  the  Reynolds-stress  tensor  is  known  as  the  “retum-to-isotropy”  term. 
It  is  part  of  the  modeled  representation  of  the  pressure/strain-rate  correlation  tensor  that  appears  in  the  exact 
differential  equation  for  the  Reynolds  stress  tensor,  r,,.  Specifically,  it  appears  as  follows. 

dr  f  7  \ 

P-^-  =  -pfi,Cla^Tg+-kSgJ  +  ...  (5) 

In  Equation  (5)  p  is  density,  t  is  time,  /?*  and  C\  are  closure  coefficients,  k  is  turbulence  kinetic  energy  and  <%  is 
the  Kronecker  delta.  Of  all  the  closure  coefficients  in  the  Stress-cy  model,  C\  is  the  least  certain.  As  argued  by 
Launder,  Reece  and  Rodi  (1975),  its  value  can  lie  within  the  following  range. 

1.4  <  Cj  <1.8  (6) 

The  Stress-co  model  uses  a  value  of  C\  =  1.8,  which  yields  optimum  results  for  a  wide  range  of  flows. 

However,  other  values  of  C\  within  the  entire  range  specified  in  Equation  (6)  yield  acceptable  results  for  the 
baseline  test  cases  that  have  been  used  to  validate  the  model.  For  each  value  of  C\  there  is  a  unique  value  of  the 
closure  coefficient  a  [see  Equation  (1)]  that  preserves  the  law  of  the  wall,  so  that  in  all  computations  we  have 
done  we  have  simultaneously  changed  both  C\  and  a. 

Disappointingly,  varying  C\  over  the  entire  range  given  in  Equation  (6)  yields  insignificant  changes.  The 
location  of  the  centerline  stagnation  point,  for  example,  changes  by  no  more  than  6%. 

4.2  k-co  Model  Revision 

Our  second  attempt  at  improving  predictive  accuracy  through  model  revisions  was  to  use  a  nonlinear  Reynolds- 
stress/strain-rate  constitutive  relation.  The  relation  we  selected  is  the  Wilcox-Rubesin  (1980)  formulation  in 
which  Zy  is  a  quadratic  function  of  the  mean-strain  rate  tensor. 

Again,  results  obtained  were  disappointing.  The  most  important  change  was  the  location  of  the  centerline 
stagnation  point.  The  use  of  the  nonlinear  stress/strain-rate  relationship  causes  the  stagnation  point  to  move 
farther  downstream,  corresponding  to  a  larger  recirculation  region.  Since  the  unmodified  k-co  model  predicts 
the  location  to  within  2%  (see  Table  1  in  Section  3),  this  is  an  undesirable  change. 

This  is  consistent  with  the  results  obtained  by  Thangam  and  Speziale  (1992)  for  flow  past  a  backward- facing 
step.  Use  of  a  nonlinear  constitutive  relationship  with  the  k-e  and  other  two-equation  models  typically  increases 
predicted  reattachment  location  by  10%. 


5.  Effects  of  Base  Bleed 


In  the  final  task  of  the  study,  we  have  analyzed  effects  of  base  bleed  on  supersonic  flow  past  a  square-based 
body.  We  consider  flow  at  two  freestream  Mach  numbers,  viz.,  Mx.  =  1.88  and  2.48  corresponding  to  the 
experiments  of  Bowman  and  Clayden  (1967).  For  both  Mach  numbers,  the  mass-injection  parameter,  /, 
varies  from  0  to  0.02,  where 


/  = 


(7) 


The  quantity  m .  is  the  mass  flow  rate  at  the  bleed  exit,  px  is  freestream  density,  Ux  is  freestream  velocity 
and  Ab  is  the  area  of  the  base. 

Figure  5  compares  computed  and  measured  area-averaged  base  pressure,  pb.  Results  of  k-e  model 
computations  [Sahu  (1986)]  are  included  for  comparison.  As  shown,  the  k-co  results  are  very  similar  to 
those  of  the  k-e  model.  For  Mach  1.88,  both  models  predict  an  increase  in  base  pressure  that  becomes  less 
rapid  as  /  increases.  For  Mach  2.48,  both  models  predict  an  initial  rise  in  base  pressure  followed  by  a  less- 
rapid-than-observed  decrease. 


P»/p°°  Pw/Poo 


Figure  5:  Computed  and  measured  effects  of  base  bleed  on  base  pressure  for  freestream  Mach  numbers  of 
1.88  (left)  and  2.48  (right).  O  Measured,  Bowman  and  Clayden;  - k-oj  model; - k-e  model. 


6.  Summary  and  Conclusions 

Results  of  this  study  and  those  of  the  preceding  study  by  Wilcox  (2001)  show  that  the  k-co  model  yields 
predictions  for  projectile  base  flows  that  are  much  closer  to  measured  properties  than  the  k-s  model.  The 
most  significant  discrepancy  between  computed  and  measured  flow  properties  is  that,  like  the  k-s  model, 
computed  base  pressure  varies  with  radial  distance  from  the  centerline.  The  area-averaged  value  predicted 
by  the  k-co  model  is  nevertheless  quite  close  to  corresponding  measurements  for  a  wide  range  of  Mach  and 
Reynolds  numbers. 

As  a  disappointing  outcome  of  the  research,  results  obtained  indicate  that,  without  additional  model 
development,  the  Stress-cn  model  offers  only  marginal  improvement  in  predictive  accuracy  over  the  k-co 
model.  In  some  cases,  the  k-co  model  predicts  flow  properties  closer  to  measurement  than  the  Stress-co 
model. 
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